Última actualización: 20 mayo 2020





Generalidades

En esta entrega se hace una verificación complementaria en función de:

  • Producción de sedimentos: Se muestra la producción anual de sedimentos para varias tipos de uso del suelo y pendientes.
  • Cantidad de agua: Se incluyen estaciones de caudal que fueron descartadas en el proceso de calibración/validación.

Sedimentos

La siguiente figura muestra un resumen de la producción anual de sedimentos ordenado según el uso de la tierra. Los usos EUCA y MONT tienen una producción de sedimentos bastante inferior AGRP, AGRC, LECH y GRAS. Adicionalmente AGRP, LECH, y GRAS tienen un comportamiento bastante similar. AGRC está ligeramente por encima y su variabilidad practicamente duplica la variabilidad del grupo AGRP, LECH y GRAS.

A continuación se muestra la producción de sedimentos en funcion de las 6 clases de pendientes modeladas en SWAT:

## [1] "S1(0-3.0)"     "S2(3.0-6.0)"   "S3(6.0-9.0)"   "S4(9.0-12.0)" 
## [5] "S5(12.0-15.0)" "S6(15.0-9999)"

Cantidad de agua

La verificación de cantidad de agua se hace a escala mensual en el período 2006-2015 (2005 de calentamiento). Esta verificacion muetra a grandes trazos que en los primeros años de simulacion (2006-2010) predomina la subestimación. Luego a partir del año 2010 se sobreestima el caudal.

A continuación leemos los resultados del modelo ejecutado en ese período así como los valores observados. Nótese que la simulación de caudal tiene de 100 ejecuciones del modelo, los estadísticos se dan en funcion de la mediana de las 100 simulaciones.

flow_m23 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/resultados/optim20200401_22/flow_m23.RDS")
flow_m30 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/resultados/optim20200401_22/flow_m30.RDS")
flow_m36 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/resultados/optim20200401_22/flow_m36.RDS")
flow_m38 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/resultados/optim20200401_22/flow_m38.RDS")
flow_m41 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/resultados/optim20200401_22/flow_m41.RDS")

flow_m23 = lapply(flow_m23[,-1],function(x){convertFlow(zoo(x, flow_m23$date),
                                                        from="cumecs", to="mm", area.km2 = 687,
                                                        timestep.default = "months")}) 
flow_m30 = lapply(flow_m30[,-1],function(x){convertFlow(zoo(x, flow_m30$date),
                                                        from="cumecs", to="mm", area.km2 = 1092,
                                                        timestep.default = "months")}) 
flow_m36 = lapply(flow_m36[,-1],function(x){convertFlow(zoo(x, flow_m36$date),
                                                        from="cumecs", to="mm", area.km2 = 2744,
                                                        timestep.default = "months")})
flow_m38 = lapply(flow_m38[,-1],function(x){convertFlow(zoo(x, flow_m38$date),
                                                        from="cumecs", to="mm", area.km2 = 3159,
                                                        timestep.default = "months")}) 
flow_m41 = lapply(flow_m41[,-1],function(x){convertFlow(zoo(x, flow_m41$date),
                                                        from="cumecs", to="mm", area.km2 = 4896,
                                                        timestep.default = "months")})

Qsim = cbind(Qsim_23 = apply(as.data.frame(flow_m23),1,median),
             Qsim_30 = apply(as.data.frame(flow_m30),1,median),
             Qsim_36 = apply(as.data.frame(flow_m36),1,median),
             Qsim_38 = apply(as.data.frame(flow_m38),1,median),
             Qsim_41 = apply(as.data.frame(flow_m41),1,median))


Qobs <- readRDS("~/R_projects/SWAT_calibration_ETmodis/data/runoff_mensual_2006_2015.RDS")

qo = cbind(Qobs, Qsim)

La siguiente tabla muestra los estadísticos KGE, NSE y Pbias a escala mensual en el período 2006-2015.

WL_23 WL_30 WL_36 WL_38 WL_41
KGE 0.25 0.32 0.74 0.37 0.70
NSE -0.26 0.34 0.73 0.33 0.79
PBIAS 47.20 -16.10 21.60 -46.50 26.90

 




A work by Rafael Navas